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, Q , Abstract 

^ ' Non Abelian vortices of a SU{2) Chern-Simons-Higgs theory in 2 + 1 dimensions are constructed numer- 

ically. They represent natural counterparts of the U(l) solutions considered by Hong, Kim and Pac, and, 
' by Jackiw and Weinberg. The Abelian embeddings are identified, for all values of the Higgs selfinteraction 

strength v, resulting in both attractive and repulsive phases. A detailed analysis of the properties of the 
solutions reveals the existence of a number of unexpected features. For a certain range of the parameter v, 
it is shown that the non Abelian vortices have lower energy than their topologically stable Abelian coun- 
^ . terparts, resulting in an effective energy lower bound on the SU{2) configurations. The angular momentum 

' of these vortices is analysed and it is found that unlike the Abelian ones, whose angular momentum and 

(~ i , energy are unrelated, there is a nontrivial mass-spin relation of the non Abelian vortices. 
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Selfdual Abelian Chern-Simons-Higgs vortices were constructed by Hong, Kim and Pac p!], and, by Jackiw 
and Weinberg [2^ describing anyonic solitons in 2 -I- 1 dimensions. These excited considerable interest because 
of their relevance to high Tc superconductivity and resulted in the development of a large body of literature, 
including also the case when the Maxwell term was present. For a complete review we refer to ^ (see also the 



, recent work d]). 
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The study of these Abelian Chern-Simons-Higgs (CS-H) vortices was motivated by the discovery in [5], of 
topologically massive (non Abelian) Yang-Mills (YM) theories augmented by a Chern-Simons (CS) term. Thus 
it is that even before the discovery of the former [Illl], their non Abelian versions were considered. Early work 
employing a non Abelian Yang-Mills-CS-H (YM-CS-H) model featuring a pair of Higgs fields in the adjoint 
representations of SU{2) was carried out in PS| PH, and in the adjoint representations of SU{N) in [S]. Non 
Abelian CS-H vortices with Higgs field in the fundamental representation of SU{2) were also discussed in [5]. 
Recently solutions in a supersymmetric JV = 2 non Abelian CS-Higgs (CS-H) theory were considered in [10], 
and more recently, vortex solutions to a U{N) YM-CS-H system with adjoint representation Higgs, describing 
the bosonic sector oi a Af — 2 supersymmetric model, were given in llj where the moduli space approach of 
[12j was employed to study the vortex dynamics at low energy. 

However, it appears that in the existing literature on this subject no numerical constructions of non Abelian 
CS-H vortices are presented. Given the absence of analytic solutions in closed form, this means that to date no 
concrete constructions are given. In the present work, we consider the simplest possible SU (2) extension of the 
2-1-1 dimensional Abelian CS-H model of [TJ[2], with adjoint representation Higgs. The Abelian embedding of 
these SU{2) vortices in the selfdual limit are none other than the selfdual vortices of [Tl[5]. 

While the physical interest of such vortices may be in condensed matter physics, here we shall put the 
main emphasis on the numerical study of these solutions, and not on their physical application. Therefore, 
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our approach is just the opposite and therefore complementary to the one generally adopted in the existing 
literature. 

Our study focuses on the full non Abelian solutions with special attention given to the relation of these 
with their Abelian counterparts. The reason is that there is no Bogomol'nyi type topological lower bound 
for the simple non Abelian system we consider, in contrast to the Abelian subsystem for which there exists 
the (magnetic) vortex number, which is the only topological charge, resulting in the topological stability of 
the Abelian vortices. It is then pertinent to examine whether the energies of the corresponding non Abelian 
vortices, namely those characterised by the same values of the physical coupling constants in the model, are 
larger or smaller than those of the Abelian embeddings. This would have a bearing on the question of the 
stability of the non Abelian vortices. To this end, we have displayed the energies of the vorticity— n solutions 
versus the Higgs selfinteraction parameter. This is an important question since these non Abelian solutions do 
not saturate a Bogomol'nyi bound, unlike for example in the presence of massive fundamental scalars [TT] . 

Another interesting question addressed is that of the value of the angular momentum of the non Abelian 
vortices, which is the only global quantity exclusively characterising these, relative to the known values for their 
Abelian counterparts [H |2] . 

Finally, as a byproduct of the present work we have constructed the non-selfdual version of the Abelian 
solutions by departing from the Bogomol'nyi limit away from the critical value of the Higgs selfinteraction 
strength. The resulting solutions exhibit the same properties of mutual attraction and repulsion that are seen 
in the usual Abelian Higgs model itself [T3] . It appears that to date this particular result has not appeared in the 
literature. After constructing these Abelian vortices, we proceed to construct the fully non Abelian solutions, 
both for the value of the Higgs selfinteraction parameter for which the Abelian embedding is selfdual, and, for 
other values of this parameter. 

The paper is structured as follows: in the next section we present the model, impose rotational symmetry 
and discuss the residual one-dimensional system. Then in section 3 we carry out the numerical constructions 
and summarise our results in section 4. 



2 The model 

2.1 The action and the general Ansatz 

The CS-Higgs model on a 2 + 1 dimensional Minkowski spacetime is described by the following Lagrangian 

C = ncs+Ti-D^<^>D^<i>^V{m,7^) (1) 

in which the CS density is 

^^cs = f e"^-- Tr Ap (^F^, - ^ A,A,^ (2) 

and the symmetry breaking Higgs selfinteraction potential is that employed in [TJ [2] 

(4A)2Tr$2(7/2 + $2)2^ (^3^ 

The dimensions of the various constants appearing above are, [77] = L~^, [A] = L and [k] = L^^ , and the index 
^ = 0, i, with i = 1,2. The Lagrangian ([l]) usually enters the more complicated models as the basic building 
block (see e.g. [TT] and the references therein). Therefore one can expect the basic features of its solutions to 
be generic. 

The static Hamiltonian of the Lagrangian ([1]) is 

Hstat = i [Tr(i^o<I'' + A<i>')+m*U)] [Tr([Ao,$]2 + A<i>')+T^(|$U)] . (4) 

We take the static "spherically" symmetric 5*0(4) YM field in 3 spacetime {i.e. 2 Euclidean) dimensions, in 
one or other chiral representation of SO±{A)^ such that our Ansatz is expressed in terms of the representation 
matrices 

sS^=-^^^)[7o,70] , a,/? =1,2,3,4, (5) 
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7q = (7i,7M)j with the index M = 3,4, being the gamma matrices in 4 dimensions and 75, the chiral matrix. 
Our rotationaUy symmetric Ansatz for the Higgs field $ and the YM connection = {Aq, Ai) is 



Ao 
A, 



(6) 



{ex)i {en)j + (eAr)^'^ Xi rij 



{sx)i 



.(±) 

^34 



in which the functions {S,^' = (x*^? X^) = X and (A*^, A^.) = Ar parametrise the YM connection in terms 



of three sets of isotriplets x ^^^d Ar, the isotriplet = <j> parametrising the Higgs field. All four 

isotriplets depend only on the 2 dimensional spacelike radial variable r, e is the two dimensional Levi-Civita 
symbol, and rii = (cos mp, sin rup) is the unit vector encoded with the winding (vortex) number > 1 in the 
{xi,X2) plane (with = xf +a:|). 

Having stated the Ansatz ([6]) in terms of the gamma matrices in 4 dimensions, we adopt henceforth the 

= (0^^ cj)^), and Ar = {Af, A^), with M = 1, 2 now. 



simpler labeling ^ = x = ix^ iX^ 



2.2 The residual system and a consistent truncation 

The parametrisation used in ^ results in 50(3) gauge covariant expressions for the YM curvature F^^ = 
{FijjFio) and the components of the covariant derivative of the Higgs field = {Di^, Dq^), expressed 

exclusively in terms of the covariant derivatives of the three triplets ^ = {^^^,£,^), X = (x*^:X^) a-^^d cf) = 
{4>^ , 0^), these covariant derivatives in the residual one dimensional space being defined as 



DrC = drC + £ A\ r 



= drXl' + e'^'^' Al X' 



ibc 



(7) 



That the residual one dimensional system of fields resulting from the imposition of this symmetry is entirely 
expressed in terms of the 5*0(3) covariant quantities ([7|) is a consequence of the consistency of Ansatz ([6]), which 
has been verified explicitly. 

The Euler-Lagrange equations arising from the variation of the gauge field are 



-e,jF,j + [$,[Ao,$]] 



(8) 



0, 

Ke^jFjo-[<^,D,<P] = 0, 
the first of which being the Gauss Law equation. The Higgs equation is 

A A$ - [Ao, [Ao, $]] - Xrf $ (772 + $2) (^2 + 3 $2^ ^ ^ (g) 
With the notation ([7|), the gauge field equations ([5]) reduce to the following set of ordinary differential equations: 



2r 
Kr 



■DrX = - 



101' X- (0.x) 



together with the constraint equation 



The Higgs equation reduces to 



1 



X Dr 



2r 



C X X- 



Dr{rDr4>)-- |er0-(0.C)^ 



+ r 



IxP - (0.x) X -\^r{v^- |0^(^;2 _ Sj^H = 0. 



(10) 



(11) 



(12) 



In (jl2p . we have rescaled the VEV 77 as r/ = i u to simplify the expression. We note that the winding number n 
in the Ansatz does not appear explicitly in the equations of motion (|10p and (|12p , nor in the contraint equation 
(|lip . It will appear only in the boundary value of the function ^'^ at the origin. 
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That (jlip is a constraint equation is easily verified by acting on it with and employing the other three 
field equations. 

In what follows we shall set the triplet of functions Ar = 0. That this is a consistent truncation is obvious 
since there is no curvature in one dimension. 

Substituting ([6]) in ([1]) yields the following reduced one dimensional Lagrangian 



L 



1, 



|2\2 



(13) 



The equations of motion (|10p and (|12p . with all covariant derivatives replaced by the ordinary derivatives 

Drc/) — > ^ = (/if, etc., follow from the variation of (jlSp with respect to the three triplets x ■> C ^-nd (f> respectively. 
It turns out that the Ansatz ([6]) can be consistently truncated further by setting 

f=(cfc^^a) , x^{dk^',b) , $^{vhk^' ,vg) , (14) 

where is a constant unit length 2— vector 0. We are thus left with only six radial functions a{r), h{r), c(r), 



(15) 



d{r), h{r) and g{r) resulting in the truncated version of the reduced Lagrangian (fT3l 
-^trunc = — — [{a br — b Qr) + {c dr — d Cp) + n br] + r (bh — dg)'^ 



r[hl + gl) + hah- cgf + \^ r (h' + g')[l - {h' + g^)f 



leading to the static energy density functional resulting from (|T5 



1 2 
- V 

2 



T {bh ~ dgf + r{hl + gl) + -{ah- cgf + r {h' + g^)[l - {h^ + g')f 

r 



(16) 



2.3 The Abelian case 

It is natural at this stage to isolate the Abelian embedding of this system, which will play an essential role in 
the construction of the non Abelian solutions. The Abelian embedding results from the truncation 

c = d = g = , 

for which the constraint ([TT|) is identically satisfied. 

The remarkable feature of this Abelian case is that the Gauss Law equation is an algebraic equation enabling 
the elimination of the electric component Aq of the Maxwell field = {Ai, Aq). This results in the reduction of 
the static Hamiltonian becoming identical with that of the usual Abelian Higgs model, subject to an additional 
proviso restricting the symmetry breaking Higgs selfinteraction potential to be of the form Q, namely the 
natural case chosen in [U [2] ■ 

The Gauss Law equation now reduces to 



2w2/i2 



leading to the reduced one dimensional Lagrangian 

Luf^i) = ^[bor- {a + l)br]+v'^rb'^ - 



(17) 



(18) 



The Gauss Law constraint p?)) . can now also be derived from the variation of p8)) with respect to 6(r). 

Next, using integration by parts, we replace the term (a + 1) br in ()18p by bar, whence ()18p can be expressed 

as 



-c/(i) 



rhi 



stat 



(19) 



^ Using diflferont constant unit length 2— vectors in ^, x a^nd </' does not lead to a consistent truncation. 
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To establish the topological lower bound of i?stat defined by , we consider the two inequalities 



a h 



^2 2y/¥h 



> 0, 



> 0, 



which lead to the final inequality 



^fstat > ±v^ [a{h^ - l)r + n\ar{h^ ~ 1)] 



(20) 



Now for this inequality to present a topological lower bound on the energy, the right hand side of (|20p must be 
a total derivative. This is only possible if we choose the constants subject to 



kA = 1 



(21) 



Saturating pO)) with (|2ip yields precisely the Bogoniol'nyi equations satisfied by the selfdual Chern-Simons 
vortices of [Tl[2], which is relevant here because our numerical analysis will depart from these vortices, first to 
the non-selfdual case analogous to the corresponding Abelian Higgs vortices [13] , and finally to the non Abelian 
vortices. 



2.4 The global charges 

On the question of the magnetic and electric fiuxes, the situation is as follows. For the generic Abelian case, 
the operative equations are the Maxwell equations, and of these the Gauss Law equation is [2] 

V ■ E - kB = p (22) 

where Ei = Fio and B — F12 ~ ^SijFij. The density p is the 0— th component of the U{1) current expressed in 
terms of the complex Higgs field such that its volume integral over cPx is the electric fiux. The volume integral 
of the divergence term V • -E, after converting it to a 'surface' (line) integral, vanishes. This can be seen readily 
from the results of the foregoing asymptotic analysis. As a result of the vanising of V • i? in (1^^ for our model, 
it is clear that the magnetic fiux J B<Px is inversely related to the electric flux. 

The only topological charge in this system is the magnetic flux, and hence also the electric flux. These are 
defined in the context of the Abelian subsystem of the SU{2) model at hand, and these respective charges are 
the only global quantities pertaining to the solutions studied. The nonabelianness therefore has no infiuence in 
this sector. There is however the angular momentum, or the spin, of these solutions, which presents another 
global quantity characterising our solutions. Unlike the magnetic and electric charges however, which are not 
infiuenced by the nonabelianess, the angular dependence is indeed dependent on the gauge group. 

The angular momentum density is 

J^T^o^ {x e), T,o = {x £)» r Tr A$ Dq^ , (23) 
which, for the fields subjected to the Ansatz ©, and further truncated according to (HH), readily yields, 

J - ^"'(^xx)-(0xe) 

= ^v^iah-cg)ibh~dg). (24) 
Remarkably enough, the total angular momentum, which is given by the integral 

J = 27r / Jrdr, (25) 
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can be expressed as a difference of two boundary integral^ (thus in this one dimensional case it reduces to the 
integral of a total derivative). It follows from the field equations ([H]), that the integral ([^5]) can be written as 

J=-«| ^.Dr^dr=-n^{\a^)\'-m\') ^^lin'-pD- (26) 

Note that pi appearing in (I26p is a asymptotic parameter appearing in ()30p below, and that the Abelian 
embedding solution is consistent with the value pi = of this parameter. Thus in this limit, the expression 
([26|) coincides with the angular momentum given in [Tl[2], i.e. J = KTrrr^/A, which holds for both selfdual and 
non-selfdual solutions. It is clear that the angular momentum of the non Abelian vortices differs essentially 
from that of its Abelian counterparts and that it can only be evaluated numerically, which will be described in 
the following section. 

We further define the energy of the solutions as the integral of (|16p. namely, 

E = 2tt [ rTQodr = 2TT [ Htmncdr . (27) 



3 Numerical results 
3.1 General features 

Although an analytic or approximate solution appears to be intractable here we present arguments for the 
existence of nontrivial solutions of the field equations (fT0| - ()12|) . 

For the full non Abelian system, the asymptotic expansion of the solutions near the origin is 

b{r) = 6o + M!^r2" + 0(r2"+2), 

K 

c{r) = -j^Aho + ^,)r-^' + 0{r-+% 
(n + 2) 2/i^u^ 

d{r) = d„r" + 0(r"+2), (28) 
h{r) = h,y' + 0(r"+2), 

All higher order terms in this expansion are fixed by the coefficient bo^d„,hn. Thus, at the origin one uses the 
following set of boundary conditions 



a|r=o = -n, drb\r=n = 0, c|r=o = rf|r=0 = ^|r=0 = , drg\r=o = . (29) 

At infinity, the finite energy requirements impose 

5 = cos a, h — sina, a = piCosa, c = pisina, d — p2sma, b = p2Cosa, (30) 

where a,pi,p2 are arbitrary constants fixed by numerics. Physically, pi and p2 corresponds to the asymptotic 
amplitudes of the effective scalars x and ^, respectively, while a somehow characterises the angle of the Higgs 
field with respect to the Abelian solution, since a = tt/2 in that limit. 

^ A similar property of the total angular momentum has been noticed in 3 + 1 dimensions for various models with gauged fields, 
see e.g. 1141 and the references there. However, in contrast to the model here, for all known d = 4 cases, the contribution to J of 
the inner boundary term vanishes. 
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Let us now concentrate on the numerical resolution of equations (fTU|) and (H^), together with the constraint 
given by equation pip . In order to do so, we employ a collocation method for boundary- value ordinary differen- 
tial equation^, equipped with an adaptive mesh selection procedure |15]. Typical mesh sizes include lO-^ — lO** 
points. The solutions have a relative accuracy of 10~^. 

It is worth noticing that equations p0 |) - p2|) may be rescaled by 

0^ , d , r ^ —r , \ - , (31) 

such that dependence on k and v disappears, remaining just a dependence on i^, which encodes the Higgs self- 
coupling parameter. For that reason, without loss of generality, we will assume in what follows that k = 2v^ and 
V — 1 and A will be rewritten as i//2. With this convention, the Abelian solutions approach the selfdual limit 
for v = I. For numerical reasons we further introduce a compactified radial coordinate defined by f = r/(l + 7'). 

After a detailed analysis of the equations, one finds that for a fixed integer n and a non-vanishing real 
the regular solutions to pO|l - (fT2|) depend on just one numerical free parameter. We have chosen it to be p2 so 
the two remaining constants in ()30|1 . namely, pi and a, are not free but are given by numerics. For that reason, 
in our numerical scheme the boundary conditions at infinity were chosen to be 

d-p2h^0 , bh-dg = , g'^ + = 1 . (32) 



To summarize, the numerical solutions are constructed by using the boundary conditions (|29[) . (|32p with the 
following input parameters: the winding number n, the Higgs self-coupling constant v and the asymptotic value 
P2 of the electric non Abelian potentiao 

Our procedure to generate non- Abelian solutions in the {n, i^, ^2} parameter space was as follows: for fixed 
integer n, one starts from the corresponding selfdual Abelian solution {h' — 1 and p2 — 0), which corresponds 
to the solution in p" , pi; moving v from 1 while keeping P2 — one generates non-selfdual Abelian solutions; 
moving p2 from zerqj while keeping 1^=1 one generates non- Abelian solutions with i/ = 1; finally, the general 
solutions are found when moving both v and p2. Also, nontrivial solutions are likely to exist for any value of 
the winding number n; however, in practice, numerics becomes more involved with increasing n. Similarly, we 
have found non- Abelian solutions for arbitrarily large values of v. 

In Figure 1 we show the functions a, 6, c, d, g, and h for a typical non-Abelian solution (with v = 2.0 
and p2 = 0.7) for n = 1. The remarkable non-Abelian nature of these solutions is inferred from the significant 
deviation from zero of the functions c, d, and g. 



3.2 The energy of the solutions 

Let us analyze the energy of the solutions. In Figure 2 we show the energy per vortex number of the Abelian 
solutions {p2 = 0) versus for n = 1, 2, 3. In the selfdual fimit (1/ = 1) the curves coincide. Below v = 1 
the energy per vortex number decreases with n for fixed v whilst above v = 1 it increases monotonically with 
n. This can be interpreted by saying that force between vortices is attractive below the selfdual limit {v = 1), 
while it is repulsive above v — I. We also include for comparison the energy per vortex number of non-Abelian 
solutions with p2 — 0.5. In that case there is no fixed selfdual limit where all the curves merge but curves cross 
in pairs at several values of v. 

The effect of non-Abelianity is exhibited in Figures 3. In Figure 3a the energy per vortex number is presented 
as a function of the parameter p2 for several values of i' for n = 1, 2, 3. Starting from the corresponding Abelian 

^Some of the solutions were also constructed by using a standard Rungc-Kutta ordinary differential equation solver. In this 
approach we evaluate the initial conditions at r = 10~^ for global tolerance 10~^^, adjusting for shooting parameters h„,bo and 
integrating towards r — > 00. We have noticed a very good agreement between the results obtained with these two different methods. 
The accuracy of the solutions was also monitored by computing a virial relation satisfied by the system l|15|l . 

*It is interesting to remark that these are also the usual input parameters for the dyonic Yang-Mills-Higgs (YMH) solutions in 
3-1-1 dimensions, see e.g. [16[ . 

^As we will show later for large values of u there are non-Abelian solutions with p2 = in addition to the Abelian ones. 
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n=l, v=2.0, P2=0.7 



J3 




0.0 0.2 0.4 0.6 0.8 1.0 
r/(l+r) 

Figure 1: A typical non-Abelian solution with v = 2.0 and p2 = 0.7 is shown for n = 1. 

P2=0.0, 0.5 




Figure 2: Energy per vortex number E/n versus the Higgs self-coupling constant v for Abelian {p2 = 0.0) and 
non-Abelian (p2 = 0.5) solutions with n = 1, 2, 3 . In the Abelian case, the curves meet for v = \, the selfdual 
limit. 



solution (p2 = 0) we mov^p2 for fixed n and v. We observe numerically that solutions exist only in the rang^ 

\P2\ <l. 

When the limit p2 = ±1 is approached the solution tends to the trivial solution 



a ~ —71 , b = —1 , c = 0, d = , g = =pl , h = Q 



(33) 



This limit requires some explanation. When \p2\ 1 the sequence of solutions tends to (j33|) in a pointwise 
way. In fact, as we can see from Figure 3a and subsequent ones, the energy for the limiting solution seems to 
depend on v and it is obviously nonvanishing. A naive computation of the energy using (j33|l however yields 



® Recall that the energy is an even function of p2 and thus similar solutions exist when p2 — > —P2 ■ 

^Physically, that means the amplitude of the electric potential at infinity is always smaller than the asymptotic value of the 
Higgs field, a feature present also in 3 -I- 1 dimensional gauged Higgs models, see e.g. [16) . 
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zero. The explanation of this apparent contradiction may be understood by analysing Figure 4. There the 
energy density (e = 27riJtrunc) of a sequence of solutions with \p2\ ^1 is presented. We observe that the energy 
density spreads and tends to zero as |p2| ^ 1 but its integral remains finite. Moreover, if one concentrates on 
functions {a, h, c, d, K\ for the sequence, one observes that for any finite value of r the sequence tends to the 
corresponding limiting value although that is not necessarily true for the value at infinity. 

For small values of v we find only one solution for each value of p2. However, for large values of v we find 
several solutions for the same value of p2 (the larger n is, the smaller v needs to be) with different value of the 
energy, in general. In Figure 3a that is clearly seen for n = 3, = 2 where in the range 0.63 < \p2\ < 0.75 three 
different non-Abelian solutions coexist for each value of p2- 

Figures 3b and 3c show the energy per vortex number versus pi and a, respectively. It is clearly seen that 
the curves look much more complicated when using pi or a as free parameter. Moreover, contrary to what 
happens when using p2 as free parameter, we do not observe any a priori bound for pi and a. 




Figure 3a: Energy per vortex number E/n versus the magnitude p2 of the electric potential at infinity for 
solutions with rt = 1, 2, 3 and several values of i^. For ly = 2,n — 5, one notices the existence of several solutions 
with the same p2- 



3.3 The angular momentum 

The angular momemtum of the solutions of Figure 3a is exhibited in Figure 5. Showing independence of v for 
the Abelian case {p2 — 0) we observe that the angular momenta of solutions with the same value of p2 but 
different values of v are different in general. As the masses of these solutions also change with changing p2 it is 
useful to consider the behaviour of the angular momentum as a function of the energy. 

That is shown in Figure 6. For small values of v, the angular momentum J is an increasing function 
of the energy E. In fact, we observe that J = solutions seem not to exist, except in the limit when the 
energy also vanishes. As v is increased, the curves develop a kink but the angular momentum still remains 
to be an increasing function of the energy. However, as v is enlarged more, one may find regions where the 
angular momentum becomes a decreasing function of the energy. Finally, for very large values of v the angular 
momentum monotonically decreases with increasing energy. This strange effect has been reported previously in 
other theories (for instance, in d = 3 + 1 Einstein-Maxwell-dilaton theory [T7], associated to counterrotation). 
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n=l 
n=2 
n=3 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 
Pi 

Figure 3b: Energy per vortex number E/n versus pi for solutions with n = 1, 2, 3 and several values of v. 
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Figure 3c: Energy per vortex number E/n versus a for solutions with n = 1, 2, 3 and several values of v. 



3.4 The issue of p2 = solutions 

Unexpectedly, the multiplicity of solutions in p2 we observe in Figure 3a as a function of z/, happens also at 
P2 — Q (although it cannot be seen there). Thus, the condition p2 — Q {i.e. \x\ ^ as r oo) does not 
characterize Abelian solutions H. These non-Abelian P2 = solutions exist for any value of the vorticity n. In 
Figure 7 we show this for solutions with n = 3. It is clear how the curves get more complicated as v increases 
with the appearance of non-Abelian p2 = solutions. Concerning this Figure 7 and the subsequent one Figure 9, 
one should notice that only half of the complete curves are presented for the sake of clarity. Since the energy is 
an even function of p2 the complete pictures would include also the mirror symmetric images, with respect to 
the P2 = line. 



°This is unique feature of d 
vanishing electric potential. 



2 + 1 CSH theories. In the better known d = 4 + 1 YMH case, |Ao| = at infinity implies a 
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n=l,v=2.0 



a 
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0.6 



0.3 



0.0 





P2=0.0 




P2=0.9 




P2=0.99 




n — n QQQ 

^j—Kl.yyy 


^^^^ r--j- 


P2=0.9999 



1 2 3 4 5 6 
log(l+r) 

Figure 4: Energy density e for a sequence of solutions with |p2| — *■ 1 for n — \ and v — 2.0. 




P2 

Figure 5: Angular momentum per vortex number J /n versus p2 for the same set of solutions of Figure 3a. 



There is also one further feature one should mention. As general rule, there are (almost alway^ non-Abelian 
solutions with lower energy than the connected Abelian solution for fixed values of n and v. This fact may take 
consequences on the stability of Abelian solutions. 

Let us analyse further these P2 = solutions. In Figure 8 we exhibit the energy i? of p2 = solutions for 
n = \. The picture is similar for any other value of n (although there are more kinks as n increases). Plotting 
the Abelian branch of solutions we observe there is a value oiv{v~ 75.0, 54.2, 55.6 for n = 1, 2, 3, respectively) 
at which the non-Abelian P2 — ^ branch branches off (represented by a large dot in Figure 8). Following the 
non-Abelian branch we observe it crosses the Abelian branch at another value of v. That means it is possible 
to find two different solutions with the same values of {n, v,p2, E}. Above that value we observe the existence 
of non-Abelian P2 = solutions with energy lower than the corresponding energy of their Abelian counterparts. 



'We have observed however some small regions where that is not the case. 
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Figure 6: Angular momentum per vortex number J/n versus energy per vortex number E/niov solutions with 
n = 1 and several values of v. 
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Figure 7: Energy E versus p2 for solutions with n = 3 and several values of v. 



3.5 Disconnected branches 

Even more strange situations may happen for large values of the Higgs self-coupling constant v. There we 

have found numerical evidence for the existence of disconnected non-Abclian branches. In general, all these 
non-Abelian solutions reported above are obtained by continuous variation of the free parameter in the range 
\pi\ < 1- However, by fine-tuning v and p2 it is possible to reach a region of the parameter space where the 
corresponding Abclian solution cannot be reached by just moving *-e., one can move p2 from -1 to 1 (passing 
through 0) without reaching the Abelian solution. An example of these disconnected branches is presented in 
Figure 9 for n = 3 and u = ZQ. 
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Figure 8: Energy E versus the Higgs self-coupling parameter ly for P2 — solutions with n = 1. For a range of 
the nonabelian solutions have a lower energy than their abelian counterparts. 
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Figure 9: Energy E versus p2 for n 
(starting point at p2 = 0) is shown . 



3, = 30 solutions. A branch disconnected from the Abelian solution 



4 Summary and discussion 

We have constructed non Abelian vortices in a SU (2) Chern-Simons-Higgs theory in 2 + 1 dimensions directly 
generalising the Abelian model proposed in [1] and [2]. These solutions are constructed numerically, and sit 
above the previously known [HE] Abelian embeddings. Our study is directed in three main directions: 

• To investigate the dependence of the energies of the non Abelian vortices on the various parameters 
characterising them. Some of these parameters, denoted as pi ,p2 and a, are the asymptotic values of 
the various fields describing the model, and since they are not independent of each other - such relations 
being only seen via the numerical process - we have chosen to focus of the most convenient one, namely 
P2, corresponding to the asymptotic value of the electric potential. In particular, the known Abelian 
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embedding solutions possess the value P2 — 0. Physically, the most significant parameter is the Higgs 
selfinteraction strength v of the symmetry breaking Higgs potential. We find it very convenient and 
interesting therefore to compare the energies of the non Abelian vortices with those of the Abelian ones 
versus the coupling i' in the P2 = case (see Fig. 8). We observe that for a range of values of v the 
energies of the non Abelian vortices are lower than those of their Abelian counterparts. But the latter 
being topologically stable, this indicates that the corresponding non Abelian vortices are also stabilised 
by this bound. 

The question of stability is a subtle one. Normally, when a new field is introduced in a classical system, 
the energy of the system becomes smaller than the original one. Examples of this are the Skyrme model 
after (diagonal) gauging with SU{2) [THj or with C/(l) Ollinij and the t/(l) gauging of the Goldstone 
soliton on [21 . In both these examples the gauged soliton is topologically stable against its own 
energy lower bound. On the other hand, when the system supporting a topologically stable soliton of the 
(purely magnetically) U{1) gauged soliton is augmented with an electric field such that the topological 
lower bound on the energy of the new system is still bounded only by the lower bound of the purely 
magnetic system, then it turns out that the energy of the electrically charged system is higher than the 
purely magnetic charged one. This is the situation in the case at hand. The Abelian subsystem has a 
topological lower bound which remains the only bound valid when the non Abelian degrees of freedom 
are introduced. One would therefore expect that the energy of the non Abelian vortex should be higher 
than the Abelian embedding. Hence we find it significant that for a certain range of the parameter v the 
energy of the non Abelian vortex is lower than that of the Abelian one. 

Apart from the question of stability illustrated by the branch structure of Figure 8, there is another inter- 
esting property of the dependence of the energy on the parameters fixing the numerical solution, illustrated 
by the branch structure of Figure 9. The latter describes a branch of solutions that is disconnected from 
the Abelian branch, plotted against the parameter p2. This phenomenon is strictly one that appears for 
large values of v. While the generic solutions are found to be connected to the Abelian embeddings, out of 
which they grow as their nonabelianness manifests itself, in the high v regime there appear disconnected 
branches. An explanation for this may be the fact that as i' grows, the contribution of the Higgs potential 
term must vanish, resulting in the constraint |0p ^ 1. Thus the dynamics changes from that of a Higgs 
model to one of a 0(3) sigma model. As it happens, the U{1) gauged embedding of this sigma model does 
indeed support topologically stable sohtons [H], so that the energy of the vortices of its 50(3) gauged 
extensions are also bounded from below, albeit with higher energy than the 50(2) gauged soliton [52]. 
This passage from a Higgs model to a sigma model has a noteworthy precedent, namely that of the 
sphalerons [23] of the standard model in the high Higgs coupling regime, where the limiting solutions 
describe the bi-sphalerons [24] coinciding with the (right) SU{2) gauged techniskyrmions [25] of the 0(4) 
sigma model. In the bi-sphaleron case, this is associated with the excitation of some extra-function in the 
ansatz of the matter fields. Therefore, we anticipate the possible existence of new solutions of the model 
considered in this work as well, which would be found beyond the truncation (|14p . 

• Given that the only topological charge in this model is the magnetic charge, or the vortex number per- 
taining to the Abelian subsystem, it is important to find some other global quantity that characterises the 
non Abelian vortices exclusively. This is the angular momentum of the vortex, which in addition to the 
globally defined magnetic and electric charges pertaining to the Abelian embedding and the energy, pro- 
vides a global quantity that characterises the non Abelian vortex. This quantity differs essentially in the 
Abelian and the non Abelian cases (see Eq. (|26p ) and gives a quantitative measure of the nonabelianness 
of the SU (2) vortices. 

The value of the angular momentum in the Abelian cases is independent of the coupling constant v. In 
that case it is also independent of the energy, whether or not the value of v is the critical one when the 
Bogomol'nyi bound is saturated as in [TJ[2]. In the non Abelian case by contrast, the angular momentum 
does depend on the energy. This yields a mass versus spin plot (see Figure 6) . Not unexpectedly, it turns 
out that the only spinless solutions are those with vanishing mass, i.e. trivial solutions with vanishing 
static energy. For a given value of u, namely for a given physical model, this spin-energy behaviour 
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is studied again, when varying our favoured parameter p2. We observe that for small values of i' the 
spin increases with energy, while for large v there exist regions where the spin decresases with increasing 
energy. Concerning this, we recall our comment above where it is pointed out that the Higgs model at 
hand becomes a sigma model in the limit of very large v. 

• Since the solutions we constructed pertain to the model in which the Higgs coupling constant p is not 
restricted to the critical Bogomol'nyi value, it is reasonable to inquire about the properties of the Abelian 
embedding vortices with respect to v. It turns out that there exist both attractive and repulsive phases 
of like charged vortices for non critical values of v for < VcriUcai and v > ^critical respectively, with non 
interacting vortices for v = identical- This is also not surprising and is the same as the situation is for the 
usual Abelian Higgs vortices. 

Finally we comment on the reason for our choice of the simplest SU{2) CS-H model, unlike the SU{N) model 
of [8] and the more general supersymmetry inspired models employed in [101 lllj . In our view, the present 
SO±{A) = SU±{2) CS-H model in 2 + 1 dimensions is the first member of a hierarchy of SO±{D + 2) CS-H 
model in I? + 1 dimensions. It is planned to study the D = 4 and 6 examples in the near future. Also, on general 
grounds, we expect the basic features of the model considered here to be generic for CS-H configurations with 
nonabelian matter fields, and thus to give an idea of the situation in a more general case. 

More immediately, we intend to revisit the problem of the present model augmented with the SU{2) YM 
term, to construct the corresponding non Abelian vortices numerically. Apart from its intrinsic value, such a 
numerical investigation would reveal some detailed properties of the analytic results of [26] . The latter work 
was carried out in the context of giving a rigorous proof of the result of Julia and Zee [27] which states that 
the vortices of the Abelian Higgs model in 2 -I- 1 dimensions cannot carry electric charge. These authors went 
further to extend the proof of the Julia and Zee theorem of the Abelian Higgs model, to the SU{2) non Abelian 
Higgs model in 2 -f 1 dimensions, i.e. that the electric charge vanishes also in that case. One would expect that 
in the limit of the CS term vanishing the (non Abelian) electric field vanishes. 
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